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ABSTRACT 



X 



The thick disc contains stars formed within the first Gyr of Galactic history, and 
little is known about their planetary systems. The Spitzer MIPS instrument was used 
to search 11 of the closest of these old low-metal stars for circumstellar debris, as a 
signpost that bodies at least as large as planetesimals were formed. A total of 22 thick 
disc stars has now been observed, after including archival data, but dust is not found 
in any of the systems. The data rule out a high incidence of debris among star systems 
from early in the Galaxy's formation. However, some stars of this very old population 
do host giant planets, at possibly more than the general incidence among low-metal 
Sun-like stars. As the Solar System contains gas giants but little cometary dust, the 
thick disc could host analogue systems that formed many Gyr before the Sun. 
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1 INTRODUCTION 

Debris discs are generated by collisions among comets or as- 
teroids, and so are a signpost of planetary systems, indicat- 
ing that bodies of at least planetesimal size formed. The de- 
bris absorbs starlight and re-emits it at longer wavelengths, 
and observations since 2003 with the Spitzer space telescope 
have yielded many new debris detections in the mid- to far- 
infrared. The relation betwe en giant planets and planetes- 
imal belts is still uncertain. iBrvden et al.l |2009t) find that 
the incidence and brightness of cool outer-system debris is 
not affected by the presence or absence of gas giants in the 
system, which is not surprising if the formation and evo- 
lution of giant planets at a few AU and comets at tens of 
AU are decoupled processes. More subtle effects may occur, 
such as the clearing out of comet belts when giant planets 
migrate, analagous to the Late Heavy Bombardment (LHB) 
era in Solar System history. This could result in brief dusty 
episodes as more planetesimals collide, and subsequent lower 
dust levels due to the ejection of many comets from the 
system. However, such LHB events seem to be atypical, as 
Spitzer has commonly found old star systems that are still 
very dusty (jBooth et al.ll2009l ). 
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iGreaves et al.l l|2007l ) argued that the incidence of giant 
planet and debris systems, and the weak overlap between 
them, could be explained if the mass of refractory elements 
at the planet-formation stage controlled the outcome. Young 
circumstellar discs with a large 'solid' mass could quickly 
form planet cores that go on to accrete gas atmospheres, 
while less substantial discs only form planetesimals by the 
time the gas dis perses and so e volve into debris systems. This 
model (see also lWvatt et all f2007])) addresses the differing 
metallicity dependences observed - metal-rich discs tend to 
have a higher solid mass and so a greater probability for 
building a planet core, while lower solid reservoirs are needed 
to bui ld comets and so debris has little dependence on metal- 
licity (iTrilling et al.ll2008l: IGreaves et al.ll2007l ; IBrvden et all 
120061 ; iBeichman et al.ll200rj| . 

Here we investigate the formation of planetary systems 
by observing some of the oldest Galactic stars with Spitzer. 
A good test of models is to observe old, metal-poor stars. If 
these proceeded only as far as planetesimal growth, a high 
incidence of debris might be seen, but perhaps modified by 
the longer time available for collisional grinding and dust 
removal com pared to mid-Galactic-age stars like the Sun 
l|Wvattll200Sl ) . The sample of stars studied here is identified 
as nearby members of the thick dis c - stars that formed 
within about the first Galactic Gyr l|Chiappini et alj|200ll ) 
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Table 1. Properties of the thick disc stars. The first two sets 
of stars form the volume-limited sample and the third set in- 
cludes those outside the complete volume or more recently iden- 
tified as thick disc objects (see text). Metallicities are for Fe ex- 
cept HD 19034 and 101259 which are combined metals. The last 
columns are the observed date and number of exposures for the 
new data (see text). 



HD [Fe/H] d type observation N exp 

(pc) date (24,70 ) 



new observations 



1 QflQ/1 


n a i 

-U.4tl 


35 


G5 


27/08/07 


o on 
z,zu 


22879 


-0.84 


21 


F9V 


03/03/07 


2,5 


29587 


-0.61 


28 


G2V 


29/10/07 


2,10 


30649 


-0.50 


30 


G1V-VI 


29/10/07 


2,15 


65583 


-0.67 


17 


G8V 


07/11/06 


2,5 


68017 


-0.42 


22 


G4V 


07/11/06 


2,5 


88725 


-0.70 


36 


G1V 


29/05/07 


2,10 


111515 


-0.52 


33 


G8V 


12/07/07 


2,20 


144579 


-0.70 


14 


G8V 


13/04/07 


2,5 


218209 


-0.43 


30 


G6V 


16/08/06 


2,20 


221830 


-0.45 


32 


F9V 


21/08/07 


2,10 






other volume-limited objects 




3795 


-0.59 


29 


K0V 






63077 


-0.83 


15 


F9V 






114729 


-0.25 


35 


GOV 






114762 


-0.70 


41 


F9V 










other thick disc stars 




6582 


-0.86 


8 


G5Vb 






23439 


-1.03 


24 


K2V.K3V 






76932 


-0.86 


21 


F9V 






101259 


-0.23 


65 


G8IV 






106516 


-0.74 


23 


F5V 






157214 


-0.41 


11 


GOV 







and have orbits inclined with respect to the Galactic Plane. 
These systems are about twice as old as the Solar System. 



2 SAMPLE SELECTION 

The targets (Tab le 1) were initially derived from 
iKaratas et ail (2005), who identified 22 Galactic thick disc 
stars near the Sun. The identification uses criteria of low 
metallicity and rotation speed, with values centred around 
[Fe/H] = -0.5 and 180 km/s respectively, removing halo 
stars that have lower v rot for the same metallicity. The sam- 
ple was then reduced to the 15 stars that lie within 40 pc of 
the Sun, to make the best use of the sensitivity of the MIPS 
instrument (|Rieke et al.| [2008). Four of thes e stars had been 
observed in previous Spitzer programs jBeichman et al.l 
l200fj ; iTrilling et 3^20081 : iBrvden et alj|2009h . and the other 
11 were observed for our Spitz er program number 30339. 
Since the survey was proposed, iReddv et al.l 1|2006D identi- 
fied further thick disc stars, several of which were observed 
in vari ous Spitzer programs, as was one of the more distant 
stars of lKaratas et alj (|2005l ). HD 101259. The total number 
of stars discussed here is 22, of which 19 have data at two 
Spitzer wavelengths. 

The (logarithmic) metallicities in this sample are - 
0.23 to -1.03, considerably lower than the mean of -0.06 
with standard deviation of 0.25 for local thin disc stars 
llValenti fe Fischerll2005l ). For nine of the volume- limited ob- 



Table 2. Observed fluxes F and errors, and excesses E,R derived 
from observed fluxes compared to the photospheres. E[24] is ex- 
pressed in magnitudes while R70 is the ratio of the observed and 
photospheric fluxes, with a significance for any dust of X70 i- e - 
excess flux above the photosphere divided by its error. 



HD 


F24 
(mJy) 


F70 
(mJy) 


E24 
(mag) 


R70 


X70 


3795 


137.7 ±0.7 


19.5 ±2.7 


0.03 


1 


.3 


1 


.5 


6582 


350 ±1 




0.05 










19034 


1 q n _|_n K 


0.1 ±2.4 


-0.02 





.0 


-0 


.9 


22879 


59 ±0.4 


15.4 ±5.2 


-0.04 


2 


.4 


1 


.7 




or; Q _|_n ^ 




-0.01 










23439B 


19.9 ±0.5 




-0.05 










29587 


38.8 ±1.2 


2.2 ±3.9 


-0.02 





.5 


-0 


.5 


30649 


48.2 ±1.2 


1 ±3.6 


0.01 





.2 


-1 


.2 


63077 


227.8 ±1.1 


15.5 ±4.1 


0.00 





.6 


-2 


.3 


65583 


65 ±0.3 


12.2 ±6.5 


-0.01 


1 


.7 





.8 


68017 


64.8 ±0.3 


6.2 ±6.1 


-0.02 





.9 


-0 


.2 


76932 


133.8 ±0.7 


14.7 ±2.1 


0.03 


1 


.0 


-0 


.0 


88725 


24.3 ±0.4 


2.3 ±3.2 


-0.03 





.9 


-0 


.1 


101259 


152.3 ±0.7 


23.3 ±2.2 


0.06 


1 


.4 


2 


.6 


106516 


80.4 ±0.4 


6.2 ±6.2 


-0.04 





.7 


-0 


.1 


111515 


20.2 ±0.5 


0.9 ±2.3 


-0.02 





.4 


-0 


.6 


114729 


62.1 ±0.7 


7.9 ±2.2 


-0.02 


1 


.2 





.5 


114762 


32.4 ±0.3 


0.2 ±3.1 


-0.05 





.1 


-1 


.1 


144579 


88.6 ±0.5 


13.1 ±4.9 


-0.02 


1 


.3 





.7 


157214 


217.5 ±1 


23.9 ±4.2 


0.04 


1 


.0 


-0 


.0 


218209 


33.4 ±1.4 


0.8 ±4.2 


-0.03 





.2 


-0 


.7 


221830 


53.5 ±0.4 


9.6 ±2.7 


-0.02 


1 


.6 


1 


.1 



jects, iTakeda et ail l|2007t) have calculated the ages based 
on isochrone fitting, a method that is most accurate for 
old stars that have evolved off the main sequence. Where 
the age is well determined (uncertainties of < 30 % for 
HD 3795, 30649, 68017, 114729, 218209 and 221830), these 
stars are dated at 9.6-11.6 Gyr old. This is in very good 
agree ment with their thick disc membership (|Karatas et al.l 
120051 ). and independent age estimates are consistent with 
this conclusion. We checked for chromospheric activity, 
and although current calibra tions of age versus activity 
ijMamaiek fe Hillenbrandll2008l . e.g.) are only for roughly so- 
lar metallicity, all of our program stars are inferred to have 
low activity consistent with ages greater than 4 Gyr. 



3 OBSERVATIONS AND DATA REDUCTION 

The eleven new stars were observed between 2006 August 
and 2007 October, using MIPS at 24 and 70 urn (Ta- 
ble 2). Exposure times were 3 and 10 seconds at these wave- 
lengths respectively, multiplied by the number of repeats 
given in Table 1. The data reducti on and instrument c alibra - 
tion procedures are described in lEngelbracht et alj l|2007l ); 
I Gordon et al.l 1120071) . The initial data reduction was by aper- 
ture photometry as described in ITrilling et all (|2008). Sub- 
sequently, these observations plus the archival data for the 
other stars were re-redu ced homogeneo usly as part of a 
Spitzer legacy catalogue (|Su et alj |2010| ) and these values 
are presented. 

To identify candidate excesses, photospheric signals 
were predicted systematically. K-band magnitudes from 
2MASS were extrapolated to 24 urn based on a set of pho- 
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□ T+B, detectable excesses 
■ T+B, actual excesses 
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Figure 1. Plot of survey sensitivity at 70 urn (T a ble 2 ) 
compared to Sun-like stars observed by iTrilling et al, 
iBeichman et al] l|2006h . denoted by T+B. The bins show rela- 
tive numbers of survey stars where, if dust was present at that 
excess level, it could be detected with 3cr confidence. The actual 
distribution of T+B dust levels is also plotted (dark blue bars). 
The excesses are above the photosphere, i.e. R70-1. 
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Figur e 2. K-[24] co lour against [Fe/H] (from NSTARS, at 
http://nstars.nau.edu) and linear fit. Data are from various 
Spitzer surveys, re-reduced homogeneously. 



tospheric solar-type stars. A few stars have K-band magni- 
tudes < 4.5 that would saturate the 2MASS images, and for 
these heritage photometry was transformed into the 2MASS 
system. These fluxes, for HD 3795, 6582, 63077, 76932 and 
157214, will be discussed in the legacy papers (Su, Rieke 
& Whitelock, in prep.); the last value is rather uncertain 
with a rms error of ~ 0.06 mag. As no 24 um excesses were 
found, the excess at 70 um was then calculated relative to 
the 24 um flux. The significance of excess at 70 um includes 
5% photometric errors plus the statistical error. 

One system, HD23439, is a binary with two resolved 
components that are treated separately here. Amongst the 
rest of the systems, HD 221830 has a binary companion at 
approximately 9 arcsec separation that could add 10-15 % 
extra photospheric flux at 70um (based on the relative K- 
band brightnesses), and HD 111515 has a similarly separated 
companion of a fainter nature, but both of these contami- 
nants are neglected here. 



4 RESULTS 

4.1 Survey limits 

From Table 2, we find no dust excess towards any of the 
stars. At 24 um, the differences between observed and pho- 
tospheric fluxes vary between -0.05 and 0.06 magnitudes, so 
there are no positive outliers. The standard deviation of the 
excesses is 0.03 magnitudes and the mean is -0.01 magni- 
tudes. At 70 um, there is one candidate excess of +2. 6a, 
but the most negative value is at — 2.3cr, so there are again 
no notably positive outliers. 

About half the stars were not detected at all at 70 um. 
This is mainly attributable to larger distances of rare thick 
disc objects compared to thin disc stars. For example, the 
median targe t dista nce is 28 pc compared to 20 pc in 
ITrilling et all (|2008h . while the average 70 um noise level 
is about 5 mjy in both cases. Three of the stars (HD 29587, 
30649, 218209) observed here also had high sky backgrounds 
that can increase the noise. Figure 1 assesses the sensitiv- 
ity achieved within the thick disc sample. Survey stars are 
binned here according to the minimum excess that could be 
detected with 3er confidence. The fraction of thick disc stars 
where only excesses above three times the photosphere can 
be detected (i.e. bins from 301-350 % upwards) is rather 
large. In the general stellar surveys, the sensitivity attained 
was generally better than this (light blue bins). However, 
many excesses above 300 % were seen in those studies (dark 
blue bins), suggesting that our survet could be adequately 
deep if the thick disc stars have typical debris . The biases 
in comparing the incidence of debris for the thin and thick 
disc stars are discussed further below. 

It was noted that some stars have 70 um fluxes well be- 
low the photospheric prediction, although with large errors. 
Since the MIPS calibration i s reliable over a wide range 
of fluxes (|Gordon et al.l [2007) , we considered whether the 
stars could be under-bright for reasons related to metallicity. 
However, for a wide sample of stars observed at 24 um, there 
is a mild trend of enhanced flux at low [Fe/H], contributing 
about 0.02 magnitudes at our median metallicity of -0.6 
(Figure 2). If the 24 um fluxes are normal or slightly high, 
it would be surprising if metal-poor stars switched to being 
under-bright somewhat longward^], at 70 um. Since stars 
observed with good signal-to-noise tend not to be under- 
bright, an optimal scaling test was carried out to check the 
observed to predicted flux ratio. The scale factor is found to 
be 1.08 + 0.06, so we conclude that there is no general 70 um 
flux deficit and the faint results are by chance. Also, the dis- 
tribution of significance values X70 (Table 2) is a slightly 
broadened Gaussian (a =1.15) centred close to zero, indi- 
cating the 70 um results are not anomalous. 



4.2 Dust excesses 

To search more deeply for cool excess among the ensemble, 
frames for the 11 new 70 um observations were stacked about 
the stellar positions, taking the resulting mean per-pixel sig- 
nals. In Figure 3, the mean photosphere is now faintly visible 



1 The stars HD 65583, 68017, 144579 and 157214 have also been 
observed with Spit zer IRS and have no excesses out to 34 |a.m 
llLawler et al]|200Sft . 
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KEY - debris discs: grey, planets: orange 
lighter colours: low-metal stars 



Figure 3. Co-added 70 urn data, showing the flux range 5-95 %. 
The averaged star signal is visible at the centre. 



at centre, although background structures from individual 
frames are also present. Aperture photometry gives a mean 
flux of 4.6 ± 1.2 mjy, consistent with minimal debris given 
the mean 70 urn photospheric flux of 5.2 mjy. 

A better estimate of the net cool dust within the whole 
thick disc sample can be found from a weighted average of 
the 70 [im excesses (Table 2). We define E70 = R70-1 and 
weight the values by 1/error 2 , where the error in E70 can 
be written as E70/x70. The net excess for the 19 stars with 
70 u.m data is 0.07 ±0.34, i.e. < 40 % above the photosphere 
at the 1-sigma upper bound. This result is however domi- 
nated by five stars with low-noise data. For a stellar effective 
temperature of around 5600 K and dust peaking in bright- 
ness at 70 [im, i.e. of temperature around 50 K, an excess of 
< 40 % gives an upper limit on fractional dus t luminosity 
Ldust/Lstar of 4 x 10~ 6 |Beichman et aL|[2006T ). 

Approximately 80 % of near by Sun-like stars are als o 
dust-free at this level individually i Greaves fc Wvattll2010F ). 
Hence, if the thin and thick disc populations have similar 
debris incidence and levels, the null result here for only 19 
stars is not surprising. In more detail, the expected number 
of debris detections at 70 fim can be estimated from Fig- 
ure 1, taking into account the limiting depths of the thick 
disc data. The number of detection per 50%-wide excess bin 
is then the fraction in the general (T+B) surveys, scaled to 
19 stars, and multiplied by the fraction of thick disc targets 
where such an excess was detectable. This yields a cumula- 
tive 2.1 expected detections, with a Poisson error of ±1.5, 
so finding no debris systems is within the errors. Similarly, 
there is a 5 % incidence of 24 urn excess for Sun-like stars 
in general (|Koerner et al.|[2009T ). so finding none among 22 
stars here is c onsistent with simila r populations. Studying 
trends by age, iGaspar et all ||200gh found a 24 [im debris- 
incidence of approximately 4.1 ± 1.5 % at 3-10 Gyr. In theo- 
retical models, some slow fading is predicted, with fractional 
dust luminosity decreasing e.g. by 25 % by 11 Gyr com- 
pared to typical m id-age Galactic systems at about 5 Gyr 
l|Lohne et al.ll20Qg|) . 

Statistics are ultimately limited by the small numbers of 
thick disc stars actually passing through the Solar vicinity at 
the present time, and in our volume-limited sample we have 
already searched for de bris among the neare st 15, or 70 %, 
of objects catalogued bv lKaratas et al] |2005l ). Fainter levels 




-3.8 -3.1 -2.4 -1.7 

log M SO | id (M Jup ) 



Figure 4. Distributions of solid masses available in the discs at 
early times (see text), for thin disc stars and the general nearby 
star population (shown with dashed boxes and darker shading). 
Giant planets can form from the three right-most bins (orange) 
and debris discs from the adjacent bin (grey); mass labels refer 
to the lowest values in the bin. 



of debris or cooler dust may be accessible to the Herschel 
satellite, now observing at wavelengths of 70 to 500 |J.m. 



5 DISCUSSION 

The debris incidence in the thick disc is either low or nor- 
mal; we can rule out a particularly high incidence of ex- 
cess in these metal-po or systems. Recent theoretical work by 
iJohansen et all (|2009r) has suggested that low-metal systems 
could have late planetesimal formation as the gas disc is dis- 
persed, in which case more collisional debris might have been 
expected around the thick disc stars, but this is not seen 
here. In contrast, some thick disc stars do host giant planets 
- an outcome thought to reflect more efficient planetesimal- 
building a t early stages. In our sample HD 114729 has a 
gas giant l|Butler et all [2003). and four other good planet- 
candidates exist for thick disc stars l|Gonzalej|2009l ). The M 
sini of these planets range from 0.5 to 7 Jupiter masses and 
the semi-major axes from 0.15 to 2 AU. 

Figure 4 sho ws the predictions of our model 
|Greaves et al.|[2007l ) based on the mass of solids available at 
the planet (esimal) formation stage. The model is based on a 
threshold of solid material needed for a successful outcome 
such as building a gas giant core, and with the solid mass 
within the disc scaling by the metal fraction in common 
with the star. The M BO nds distribution is obtained by mul- 
tiplying the log-normal dis tribution of overall disc masses 
l|Andrews fc Wi lliams 2005|) by the log-normal distribution 
[Fe/H] for the population of interest, here with mean and 
standard deviation of -0.6 and 0.2 respectively. The thresh- 
olds for forming detectable debris discs and gas giants are 
around 0.02 and 0. 1 Jupiter masses in solids respectively 
|Greaves et al.f 2007). For simplicity we neglect systems with 
both planets and debris, since t he two phenomena appear 
uncorrelated l|Brvden et al.|[2009l ). and so this is a small sub- 
set of Sun-like stars. 
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Figure 4 then predicts that 15 % of thick disc stars 
should form giant planets, with a further 9 % forming comet 
belts. Around the general local population of Sun-like stars, 
the corresponding numbers are 21 % and 20 %. The differ- 
ences are modest because the availability of solid material is 
largely dictated by the total disc masses, which have a very 
broad distribution of a = 1.15 dex, compared to the stan- 
dard deviations around 0.2 dex for the metallicities. Hence, 
while the peaks of the solid-mass distributions do shift for 
differing metals (Figure 4), the upper tails associated with 
forming planet (esimals) are broad and do not have greatly 
different numbers of stars. 

These values are compatible with survey data, taking 
into account that not all systems of planets and debris have 
yet been discovered. In the general stellar population, the 
extrapolated planet total is 17-20 % (|Cumming et al.|[2009T ) 
and the debris disc count at 70 urn is aro und 14-20 % 
l|Trilring et all |200S| ; iGreaves fc Wvatdl2010l ). which agree 
well with the predictions of 21 and 20 % respectively. For the 
thick disc stars, the debris incidence is < 5 % versus the 9 % 
predicted, and the planet count is estimated at ~ 6 — 11 %, 
versus 15 % predicted. It has been noted t hat lower- metal 
discs may disperse faster |Yasui et al.ll2009h , and this could 
reduce the efficiency of growing planet (esimals) compared 
to our simple model. 

The thick disc planet frequency is poorly known, 
with our 6 % estimate derived from three planetary sys- 
tems among 52 thick disc objects that have [Fe/H] < 
-0.2 and [a/Fe] > 0.2 l|Reid et al.ll2007l ; IValenti fc Fischer! 
2005), while the 11 % is for two out of 18 systems in 
the Fe,Mg-selected, volume- limited, thick disc sample of 
iFuhrmann fc BernkopJ |2008). These planet frequencies are 
actually higher than for low-metal stars in general, such as 
the incidence of 1 % f or periods of less than three years found 
in a recent survey bv lSozzetti et all |2009r i. Conditions that 
could promote the formation of more gas giants in the old 
thick disc population compared to younger low-metal thin 
disc stars are at present unclear. 



6 CONCLUSIONS 

While the scarcity of nearby thick disc stars makes it diffi- 
cult to study the population, there are examples of disc and 
planet systems around some comparably old stars. A gas gi- 
ant is known aro und the post main s equence star V391 Peg 



of age > 10 Gyr (|Silvotti et al 



20071); the w 10 Gv r-old G8 



^yr 

V star r Ceti has a debris disc ( Greaves et al.ll2004l ) ; and the 



oldest know n star with both a planet and debris is 70 Vir 
at w 9 Gyr (|Brvden et al.ll2009h . These discoveries point to 
some of the first stars in the Galaxy having had the poten- 
tial to form planetary systems akin to those around typical 
stars seen today. If some thick disc stars have minimal cool 
debris but do host gas giants orbiting at a few AU, these 
would be analogues to the Solar System formed many Gyr 
before the Sun. 
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